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for  the  proposed  Arcadia  Lake. 
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Preface 

The  work  described  in  this  report  was  performed  by  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi,  for 
the  U.  S.  Army  Engineer  District,  Tulsa.  The  project  was  authorized  by 
Intra-Army  Order  for  Reimbursable  Services  No.  DOL  7^0015,  dated 
2 August  1971*  and  amended  31  August  1976  for  the  present  sJudy. 

This  report  describes  an  evaluation  of  the  water  quality  expected 
in  the  proposed  Arcadia  Lake  at  the  lower  pool  elevation,  305  m msl. 

This  report  is  intended  to  supplement  the  principal  report,  "Arcadia 
Lake  Water-Quality  Evaluation,"  conducted  at  pool  elevation,  311^  msl. 

The  project  was  undertaken  by  the  Environmental » Effects  Laboratory 
(EEL)  at  the  WES.  The  research  was  conducted  under  the  direct  super- 
vision of  Mr.  D.  L.  Robey,  Chief,  Ecosystem  Modeling  Branch,  and  the 
general  supervision  of  Drs.  R.  L.  Eley,  Chief,  Ecosystem^Research  and 
Simulation  Division,  and  John  Harrison,  Chief,  EEL.  Drs.  K.  W.  Thornton 
and  D.  E.  Ford  served  as  principal  investigators.  Mr.  R.  R.  Hall,  Jr., 
participated  in  the  study  and  assisted  in  the  preparation  of  this  report. 

Director  of  the  WES  during  the  preparation  and  publication  of  this 
report  was  COL  J.  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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WATER-QUALITY  EVALUATION  OF  A LOWER  POOL  ELEVATION 
FOR  PROPOSED  ARCADIA  LAKE,  OKLAHOMA 


Introduction 


1.  This  report  presents  results  of  a water-quality  evaluation  of 
the  proposed  Arcadia  Lake  at  pool  el  305  m msl  based  on  mathematical 
ecological  modeling  studies.  This  report  was  prepared  to  provide  the 
U.  S.  Army  Engineer  District,  Tulsa,  with  guidance  on  the  impact  of  a 
lower  pool  elevation  on  reservoir  water  quality. 

2.  Another  more  intensive  study,  "Arcadia  Lake  Water-Quality 
Evaluation,"'*'  was  conducted  at  pool  el  311  m msl.  This  study  supplements 
the  other  analyses.  It  is  recommended  that  the  "Arcadia  Lake  Water- 
Quality  Evaluation"  report  be  consulted  for  specific  details  concerning 
the  mathematical  modeling  approach.  For  the  remainder  of  this  report, 
the  study  conducted  at  pool  el  311  m will  be  referred  to  as  the  principal 
study. 

3.  The  project  as  authorized  includes  a multiple-purpose  reservoir 

located  on  the  Deep  Fork  River  near  Arcadia,  Oklahoma,  for  flood  control, 

municipal  and  industrial  water  supply,  and  recreation.  The  project,  as 

formulated  for  this  study,  would  have  a conservation  pool  IT  m deep  near 
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the  dam  at  el  305  m and  would  impound  approximately  3.^  x 10  m of 
water  having  a surface  area  of  736  ha.  This  would  create  a lake  extend- 
ing approximately  9 km  upstream  from  the  dam  with  a mean  depth  of  4 m. 

Simulation  Approach 


U.  The  water-quality  simulations  were  conducted  using  an  improved 
version  of  the  Water  Quality  for  River-Reservoir  Systems  (WQRRS)  reser- 
voir model.  Several  modifications  and  improvements  in  the  WQRRS  model 
have  been  made  since  it  was  first  applied  to  Arcadia  Lake  in  1975. 
Results  should  be  analyzed  in  terms  of  relative  increases  or  decreases 
in  magnitude  and  phasing  and  not  in  absolute  terms. 
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5.  Chemical  and  biological  modeling  coefficients  were  identical  to 
those  used  in  the  principal  study.  Modifications  to  the  physical  com- 
ponents of  the  model  required  changes  in  several  physical  coefficients. 
Simulations  were  conducted  using  inflows  previously  modified  to  repre- 
sent the  diversion  of  the  Northside  Sewage  Treatment  Plant  effluent  from 
the  upper  Deep  Fork  River.  Model  simulations  used  the  same  daily  mete- 
orological updates  and  monthly  streamflow,  physical.,  chemical,  and  bio- 
logical updates  for  study  years  1970  and  1973  as  were  used  in  the 
principal  study.  These  two  years  represented  below-  and  above-average 
flow  years,  respectively,  for  the  period  from  1969  to  197*+  for  which 
stream  water  quality  data  were  available.  The  chemical  updates  had 
previously  been  modified  to  consider  the  diversion  of  the  sewage 
treatment  plant. 


Water-Quality  Simulations 


6.  Table  1 lists  a series  of  simulations,  and  their  objectives, 

that  were  run  on  the  1970  and  1973  data  sets.  The  data  from  these  simu- 
lations were  used  to  investigate  the  sensitivity  of  the  model  to  various 
modeling  coefficients  (Table  l).  Since  there  is  some  error  associated 
with  each  model  coefficient  or  rate  constant  similar  to  the  errors 
associated  with  sampling  or  measuring  river  inflows  (model  updates), 
the  model's  sensitivity  to  changes  in  coefficients  and  updates  was 
investigated  to  determine  if  this  variability  would  significantly  affect 
conclusions  drawn  from  the  modeling  results.  The  results  are  discussed 
below  for  the  following  conditions:  (a)  base  case,  (b)  varying  coeffi- 

cients, (c)  varying  light  penetration,  (d)  varying  nutrient  inflow,  and 
(e)  el  308  m. 

Base  case 

7.  The  reservoir  was  predicted  to  thermally  stratify  intermit- 
tently during  the  spring  and  summer  months  but,  in  general,  was  iso- 
thermal at  the  lower  pool  elevation  (Figures  1 and  2).  During  both 
1970  and  1973,  wind  mixing  was  a dominant  factor  affecting  the  thermal 
structure  of  the  impoundment.  The  reservoir  would  begin  to  stratify 
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during  periods  of  warm,  calm  weather  hut  would  mix  during  periods  of 

windy  weather.  This  is  typical  of  Oklahoma  impoundments.  The  same 

pattern  of  intermittent  thermal  stratification  was  found  to  occur  in 

2 

Keystone  Reservoir,  Oklahoma. 

8.  The  isothermal  conditions  in  the  reservoir  resulted  in  release 
temperatures  similar  to  reservoir  temperatures.  Target  temperatures  for 
releases  were  selected  to  correspond  to  natural  river  temperatures  and 
were  obtained  from  long-term  monthly  averages.  In  general,  the  release 
temperatures  were  lower  than  target  temperatures  during  the  winter  and 
early  spring  and  higher  than  target  temperatures  during  the  summer  and 
early  fall  in  1970  (Figure  3).  Meteorological  conditions  during  1973 
resulted  in  reservoir  temperatures  increasing  more  slowly  in  the  spring 
but  staying  warmer  in  the  fall.  This  was  reflected  in  release  tempera- 
tures that  lagged  target  temperatures  during  the  spring  and  early  summer 
but  were  greater  than  target  temperatures  during  the  fall  and  early 
winter  (Figure  1+).  Because  of  the  generally  isothermal  conditions,  the 
selective  withdrawal  capability  for  meeting  release  temperature  objec- 
tives was  not  fully  realized.  However,  release  temperatures  always  were 
within  about  7°C  of  target  temperatures. 

9-  The  dissolved  oxygen  (DO)  budget  was  influenced  by  periods 
of  wind  mixing.  In  general,  DO  was  greater  than  5 mg/ SL  throughout  the 
pool  during  most  of  1970  (Figure  5).  During  periods  of  intermittent 
thermal  stratification,  hypolimnetic  DO  concentrations  decreased  below 
5 mg /l  and,  in  three  instances,  reached  0 mg/ £ in  the  bottom  layers. 

The  average  duration  of  the  stratification  periods  was  approximately  six 
days.  DO  concentrations,  in  general,  were  lower  during  1973  than  during 
1970  (Figure  6).  DO  concentrations  were  also  less  than  5 mg/ l more 
often  and  for  longer  periods  of  time.  Phytoplankton  concentrations  were 
greater  in  1973,  resulting  in  a greater  DO  demand  from  respiration 
and  the  subsequent  decay  processes.  In  addition,  1973  had  greater  in- 
flows than  1970  and,  therefore,  higher  detrital  loadings  that  also  used 
DO  during  their  decay. 

10.  There  were  two  major  algal  blooms  predicted  for  1970  (Fig- 
ure 7)*  The  first  algal  bloom  consisted  primarily  of  nannophytoplankton 
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and  occurred  from  April  through  May.  Nannophytoplankton  were  considered 
to  he  less  than  63  ym  in  size  and  would  pass  through  a No.  20  plankton 
net.  The  second  bloom  was  composed  of  nearly  equal  concentrations  of 
net-  and  nannophytoplankton  and  lasted  from  early  July  through  September. 
Three  major  phytoplankton  blooms  occurred  during  1973  (Figure  8).  In 
general,  the  magnitude  of  the  algal  blooms  was  greater  during  1973  than 
in  1970  and  all  blooms  consisted  primarily  of  nannoplankton.  The  first 
bloom  began  in  March  and  continued  until  the  end  of  May.  The  second 
bloom  occurred  from  early  June  to  mid-July,  while  the  third  bloom  began 
in  early  August  and  continued  until  the  end  of  the  simulation  period. 
Varying  model  coefficients 

11.  The  values  specified  for  each  coefficient,  in  general,  repre- 
sented a mean  or  average  value  for  the  rate  process.  These  values  were 
obtained  from  the  literature,  personal  experience,  or  personal  communica- 
tion with  authorities  on  that  particular  rate  process  in  Oklahoma  im- 
poundments. Because  there  was  a variance  associated  with  each  coeffici- 
ent, a series  of  sensitivity  analyses  were  conducted  on  coefficients 
affecting  a major  compartment  of  interest  such  as  phytoplankton.  The 
sensitivity  analyses  permitted  an  assessment  of  how  water  quality  con- 
clusions might  have  changed  if  different  coefficient  values  had  been 
selected. 


12.  Sensitivity  analyses  were  conducted  on  Algae  1 light,  nitro- 
gen and  phosphorus  half-saturation  coefficients.  Algae  2 light  half- 
saturation coefficient,  and  the  self-shading  coefficient.  None  of  the 
coefficient  changes  had  a significant  effect  on  thermal  stratification. 
In  addition,  the  nitrogen  and  phosphorus  half-saturation  coefficients 
had  no  apparent  effect  on  any  of  the  simulation  results.  Changing  the 
light  half-saturation  and  the  self-shading  coefficients  had  minimal  ef- 
fect on  the  DO  concentrations  for  the  two  study  years.  During  1970, 
varying  light  half-saturation  coefficients  resulted  in  slightly  lower 
DO  concentrations  during  mid-May  but  slightly  greater  DO  concentrations 
in  early  July  (Figure  9).  Overall,  the  DO  budget  was  quite  similar. 
During  1973,  the  annual  DO  budget  was  slightly  greater  by  varying  the 
light  half-saturation  coefficient  (Figure  10).  Increasing  the 
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self-shading  coefficient  resulted  in  slightly  higher  DO  concentrations 
during  both  1970  and  1973  (Figures  11  and  12). 

13.  Varying  the  light  half-saturation  coefficient  had  a signifi- 
cant effect  on  predicted  algae  concentrations.  During  1970,  varying 
the  algae  light  half-saturation  coefficients  resulted  in  the  net  phyto- 
plankton becoming  the  dominant  community  (Figure  13).  The  spring  bloom 
was  composed  of  approximately  two-thirds  net  plankton  and  one-third 
nannoplankton,  while  the  late  summer  bloom  consisted  almost  entirely  of 
net  phytoplankton.  The  phytoplankton  also  maintained  higher  concentra- 
tions between  blooms.  Varying  the  algae  light  half-saturation  coeffi- 
cients during  1973  resulted  in  slightly  lower  blooms  during  early 
spring  and  summer  with  a significant  shift  in  dominance  within  the  plank- 
ton community  during  late  summer  (Figure  lU) . The  late  summer  bloom 
consisted  of  an  initial  peak  of  net  plankton  followed  by  a bloom  of  sim- 
ilar magnitude  of  nannoplankton.  The  magnitude  and  duration  of  this 
late  summer  bloom  was  significantly  less  following  the  change  in  the 
light  half-saturation  coefficients. 

lU.  During  1970,  the  self-shading  coefficient  resulted  in  a 
lower  spring  bloom  but  a higher  summer  bloom  than  the  base  (Figure  15). 
Light  limited  the  increase  during  the  first  bloom.  However,  during  the 
second  bloom,  zooplankton  grazing  on  the  plankton  community  was  lower, 
resulting  in  a greater  summer  bloom.  Since  the  first  bloom  was  not  as 
large  with  the  greater  self-shading,  the  zooplankton  population  did 
not  increase  in  magnitude  as  it  did  during  the  base  case  and  initially 
was  not  able  to  graze  as  much  on  the  phytoplankton  community  during  June. 
The  second  bloom,  therefore,  began  earlier  and  was  of  greater  magnitude 
at  the  higher  self-shading  value.  Solar  radiation  was  also  greater 
during  the  second  bloom  than  during  the  first  bloom.  During  1973,  the 
blooms  were  delayed  slightly  at  the  higher  self-shading  values  but  were 
similar  in  magnitude  to  the  blooms  in  the  base  case  (Figure  l6). 

Light  penetration 

15.  The  effect  of  increased  and  decreased  light  penetration  on 
water  quality  was  investigated  by  varying  the  Secchi  disk  depth.  In 
the  WQRRS  model,  the  Secchi  disk  depth  is  used  as  a relative  measure  for 
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the  depth  of  light  penetration  into  the  reservoir. 

16.  The  Secchi  disk  depth  had  no  significant  effect  on  thermal 

1 

stratification  and  minimal  effects  on  the  DO  budget.  In  general,  the 
DO  budget  was  slightly  greater  at  the  lower  Secchi  disk  value  of  1 m 
than  the  base  case  for  both  study  years  (Figures  17  and  18) . The  DO 
budgets  for  the  base  case  and  higher  Secchi  disk  values  of  3 m were 
similar  throughout  the  year  for  both  1970  and  1973  (Figures  19  and  20). 

17.  During  1970,  the  first  bloom  was  delayed  by  approximately 
one  month  at  the  lower  Secchi  disk  value  of  1 m,  while  the  second  bloom 
was  delayed  by  approximately  15  days  (Figure  21 ).  Although  the  second 
bloom  was  delayed  at  the  Secchi  disk  of  1 m,  the  magnitude  of  the  bloom 
was  slightly  greater.  The  first  bloom  was  less  at  the  Secchi  disk  of 

1 m.  During  1973,  all  algae  blooms  at  a Secchi  disk  of  1 m were  delayed 
slightly  in  starting  but  the  magnitude  of  the  blooms  was  similar  for  all 
Secchi  disk  values  (Figure  22). 

Varying  nutrient  inflow 

18.  The  simulated  update  nutrient  concentrations  (NH^-N,  NO^-N, 
and  PO^-P)  flowing  into  Arcadia  were  doubled  and  halved  to  evaluate 
their  impact  on  reservoir  water  quality.  The  algae  blooms  were  iden- 
tical to  the  base  case.  Since  light  was  limiting  the  phytoplankton 
assemblage,  doubling  the  nutrients  had  no  effect  on  the  phytoplankton. 

It  appears  that  advanced  waste  treatment  would  have  to  reduce  nutrient 
concentrations  much  lower  than  50  percent  of  their  present  values 
before  any  apparent  changes  might  occur  in  the  reservoir  phytoplankton 

/ 

assemblage. 

Elevation  308 

19.  Since  some  quality  constituents  in  the  present  study  increased 
while  others  decreased  tfhen  compared  to  the  principal  study,  simulations 
were  therefore  conducted  to  determine  if  an  intermediate  elevation  would 
have  ameliorated  these  conditions. 

20.  The  thermal  stratification  at  the  higher  pool  elevation 

was  similar  to  the  base  case  for  1970  and  1973  (Figures  23  and  2U).  The 
reservoir  was  predicted  to  intermittently  stratify  during  the  spring  and 
summer  months  but,  in  general,  was  isothermal.  Although  the  reservoir 
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elevation  was  increased  by  3 m,  wind  mixing  was  still  a dominant  factor 
affecting  the  thermal  structure  of  the  impoundment. 


21.  The  DO  concentrations  were  greater  at  the  higher  elevation 
than  the  base  case  for  both  1970  and  1973  (Figures  25  and  26).  The 
greater  volume  of  water  and  the  frequent  mixing  resulted  in  a larger 

DO  mass  available  to  satisfy  a similar  DO  demand  to  the  base  case.  The 
net  result,  therefore,  was  an  increase  in  the  overall  DO  budget. 

22.  The  spring  plankton  blooms  were  similar  at  the  lower  and 
higher  pool  elevation  during  1970,  but  the  magnitude  of  the  summer 
nannoplankton  bloom  was  greater  at  the  higher  pool  elevation  (Figure  27). 
The  plankton  blooms  were  all  diminished  at  the  higher  pool  elevation 
during  1973  (Figure  28).  The  timing  of  the  three  blooms  was  similar 

at  both  elevations. 

Comparison  with  Principal  Study 

23.  Simulations  at  el  311  m indicated  the  reservoir  would  be 
weakly  stratified  from  mid  to  late  spring  through  late  summer  or  early 
autumn.  The  intermittent  stratification  predicted  at  the  lower  pool 
el,  305  m,  was  expected  considering  the  relatively  high  average  wind 
velocities  in  Oklahoma.  Although  complete  mixing  in  the  summer  occurred 
less  frequently  in  the  principal  study,  stratification  also  was  dis- 
rupted during  extended  periods  of  high  wind  velocity. 

24.  Intermittent  rather  than  prolonged  stratification  at  the 
lower  pool  elevation  did  result  in  higher  DO  than  in  the  principal  study. 
Wind  mixing  continuously  replenished  the  DO  depleted  by  respiration  and 
decay  processes  in  the  hypolimnion.  During  brief  periods  of  stratifica- 
tion, the  hypolimnetic  DO  decreased  and,  at  times,  the  hypolimnion  did 
become  anoxic.  During  the  principal  study,  however,  prolonged  stratifi- 
cation occurred,  a greater  volume  of  the  hypolimnion  became  anoxic,  and 
anoxia  persisted  for  longer  periods  of  time. 


25.  The  first  algae  bloom  was  similar  in  magnitude  and  duration 
to  the  bloom  in  the  principal  study  for  1970,  but  the  second  bloom  was 
lower  in  magnitude  and  slightly  longer.  The  magnitude  of  the  first  and 
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third  blooms  during  1973  was  nearly  double  the  concentration  found  in 
the  principal  study.  Net  phytoplankton  contributed  more  to  the  second 
bloom  during  1970  at  the  lower  pool  elevation  than  in  the  principal 
study  but  was  very  similar  in  concentration  to  the  principal  study 
during  1973. 

26.  State  of  Oklahoma  water-quality  standards  and  Environmental 
Protection  Agency  criteria  for  fecal  coliform  bacteria  were  not  violated 
at  any  time  or  at  any  depth  during  the  1970  or  1973  simulations  at 
either  pool  elevation.  The  simulations  are  considered  most  representa- 
tive of  the  area  in  the  deepest  part  of  the  pool  near  the  dam.  Simula- 
tion results  should  not  be  interpreted  to  mean  violations  may  not  occur 
in  the  headwater  areas  or  in  isolated  coves.  Violations  during  major 
storm  events  (storms  with  2-  to  5-day  hydraulic  residence  times  in  the 
pool)  may  occur  throughout  the  impoundment.  During  most  storm  events, 
violations  may  be  expected  to  occur  in  the  upper  third  of-  the  proposed 
impoundment . 

Discussion 


27-  It  is  important  to  realize  that  mathematical  ecosystem  models 
are  in  the  early  stages  of  development  and  verification.  The  original 
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documentation  of  the  WRE  model  was  published  by  Chen  and  Orlob  in  1972. 
During  the  succeeding  five  years,  there  have  been  significant  improve- 
ments made  to  the  model,  but  neither  the  WQRRS  model  nor  any  other  mathe- 
matical water-quality  model  is  capable  of  predicting  so-called  "absolute 
values."  Preliminary  research  evaluations  indicate  that  predictive  val- 
ues for  some  parameters  at  times  may  be  in  error  by  more  than  an  order 
of  magnitude  of  measured  values  or  by  two  weeks  of  the  onset  of  a spe- 
cific event.  However,  if  the  assumptions  and  limitations  of  the  model 
are  understood  and  incorporated  in  the  interpretation  of  the  output, 
the  general  response  predicted  by  the  model  will  provide  valuable  in- 
sight into  the  water  quality  of  the  proposed  Arcadia  Lake.  Despite 
their  limitations,  existing  modeling  techniques  properly  used  in  con- 
junction with  other  appropriate  study  approaches  may  be  considered 
as  one  of  the  more  reliable  and  versatile  tools  for  making  water- 
quality  analyses. 
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28.  The  WQKRS  model  assumes  that  a reservoir  consists  of  a series 
of  well-mixed  horizontal  layers.  Thus,  while  spatial  concentration 
gradients  in  the  vertical  direction  can  he  evaluated,  spatial  gradients 
cannot  he  evaluated  in  the  horizontal  direction.  The  one-dimensional 
assumption  of  the  WQRRS  is  important  in  the  interpretation  of  modeling 
results.  The  predicted  water-quality  conditions  are  considered  most 
representative  for  the  part  of  the  pool  near  the  dam.  This  area  has  the 
greatest  volume,  the  most  representative  stratification,  and  the  most 
significant  impact  on  downstream  releases.  The  one-dimensional  assump- 
tion, however,  implies  that  the  model  cannot  predict  isolated  water- 
quality  conditions  in  the  headwaters,  coves,  or  emhayment  areas.  These 
shallow,  usually  well-mixed  areas  cannot  he  represented  hy  a one- 
dimensional model. 

29.  The  simulations  also  are  considered  representative  of  condi- 
tions six  to  eight  years  after  filling  of  the  reservoir  when  it  is  in  a 
less  transient  state.  Transient  biological  and  chemical  conditions  are 
known  to  occur  during  the  first  several  years  after  reservoir  filling. 
These  transient  conditions  were  not  considered  in  the  simulations  con- 
ducted as  part  of  this  study.  With  an  understanding  of  these  basic  as- 
sumptions, various  predicted  water-quality  conditions  can  be  discussed. 

30.  The  predicted  release  temperatures  and  thermal  stratification 
profiles  are  probably  quite  representative  of  actual  values  that  will 
occur  in  Arcadia  Lake  if  the  project  is  constructed.  The  thermal  pre- 
dictions are  based  on  physical  processes  that  are  reasonably  well 
understood.  The  impoundment  is  predicted  to  be  intermittently  strati- 
fied. Therefore,  the  temperature  structure  will  be  highly  dependent  on 
local  hydrometeorological  conditions.  The  intermittent  stratification 
would  also  make  it  difficult  to  meet  the  average  natural  stream  release 
temperature  objectives  during  the  summer  and  early  fall  months.  In 
contrast,  at  the  higher  elevation,  the  summer  stratification  was  pre- 
dicted to  be  somewhat  stable  and  the  downstream  release  objectives  could 
be  met  through  the  operation  of  the  selective  withdrawal  structure. 
However,  even  during  intermittent  stratification,  release  temperatures 
can  be  maintained  within  T°C  of  average  natural  stream  temperatures. 
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31.  The  seasonal  DO  concentrations  predicted  for  the  proposed 
Arcadia  Lake  are  dependent  on  the  thermal,  structure,  mixing,  and  the 
phasing  and  magnitude  of  algal  blooms.  While  no  prolonged  periods  of 
stratification  were  predicted  during  1970  or  1973,  a stable  period  of 
stratification  could  occur  during  a prolonged  period  of  warm,  calm 
weather.  Under  these  conditions  it  is  probable  that  DO  would  go  to  zero 
in  the  deeper  strata  of  the  impoundment.  Bottom  withdrawal,  in  com- 
parison with  surface  withdrawal,  would  be  more  effective  in  keeping 
DO  concentrations  higher  in  the  deeper  strata  by  withdrawal  of  the  low  DO 
water  and  replacement  of  it  with  higher  DO  water  from  above  or  from  in- 
flowing waters.  Reaeration  through  the  structure  and  stilling  basin 
would  result  in  80  to  95  percent  DO  saturation  of  the  release  waters 
even  during  withdrawal  of  anoxic  bottom  waters.  This  represents  physical 
reaeration,  however,  that  may  not  satisfy  the  oxygen  demand  in  the  re- 
lease waters,  resulting  in  an  oxygen  sag  farther  downstream.  Based  on 
an  analysis  of  DO  data  downstream  from  other  Oklahoma  reservoirs  when 

anoxic  bottom  waters  were  released,  the  oxygen  sag  would  not  be  expected 
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to  violate  stream  DO  objectives. 

32.  As  indicated  by  the  sensitivity  analyses  and  light  penetra- 
tion simulations,  light  would  be  an  important  factor  regulating  plankton 
growth.  The  lower  pool  elevation  with  the  smaller  volume  is  expected  to 
be  more  turbid  for  two  reasons.  First,  wind  mixing  will  have  a greater 
mixing  effect  at  the  lower  pool  elevation  than  in  the  previous  study. 

This  should  result  in  a greater  suspension  of  matter  throughout  the  water 
column  and  possible  resuspension  of  bottom  particles  during  periods  of 
high  wind  velocity. ^ Second,  the  theoretical  hydraulic  residence  time 
has  been  decreased  by  nearly  one  third  in  both  years  (3.2  to  1.2  years 
and  1.6  to  0.5  years  during  1970  and  1973,  respectively).  The  residence 
time  during  1973  is  now  predicted  to  be  less  than  1 year.  The  higher 
flushing  rates  at  the  lower  elevation  should  generally  result  in  greater 
in-pool  turbulence,  which  would  tend  to  suspend  particles  for  longer 
periods  of  time.  This  increased  turbidity  may  reduce  the  magnitude  of 
the  predicted  phytoplankton  populations  but  it  may  result  in  coliforms, 
heavy  metals,  and  pesticides  adsorbed  onto  sediment  particles  to  remain 
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in  the  water  column  for  longer  periods  of  time  before  settling  out. 

Turbidity  and  suspended  solids  concentrations  are  not  currently  simu- 
lated with  the  WQRRS  model.  However,  the  sensitivity  analyses  does 
indicate  there  would  be  a fairly  large  range  of  algal  responses  depend- 
ing on  the  light  regime.  The  magnitude  of  the  phytoplankton  blooms  is 
expected  to  be  on  the  order  of  that  occurring  in  many  Oklahoma  reservoirs 
such  as  Thunderbird  Lake,  Keystone  Reservoir,  and  Lake  Carl  Blackwell. 

33-  Doubling  or  halving  the  inflowing  nutrient  concentrations  or 
varying  the  nitrogen  or  phosphorus  half-saturation  coefficients  did  not 
influence  the  phytoplankton  because  the  inflowing  and  in-lake  nutrient 
concentrations  were  over  an  order  of  magnitude  higher  than  those  levels 
considered  necessary  to  support  phytoplankton  growth.  Since  the 
Northside  Sewage  Treatment  Plant  effluent  was  considered  to  be  diverted 
outside  the  basin,  the  nutrients  are  those  resulting  from  nonpoint 
sources  that  are  difficult  to  control  and  reduce.  The  model  simulations 
indicated  that  substantial  reductions  in  nutrient  concentrations  would 
be  required  before  any  significant  decrease  in  algae  blooms  occurred. 

3^.  Simulations  indicate  that  the  reservoir  would  be  intermit- 
tently stratified  at  the  lower  pool  elevation  of  305  m.  In  general, 

DO  concentrations  would  be  greater  than  at  pool  el  311  m.  The  aerobic 
environment  in  the  hypolimnion  should  reduce  the  level  of  iron  and 
manganese  near  the  dam  below  the  level  predicted  to  occur  in  the  princi- 
pal study.  However,  manganese  concentrations  probably  would  exceed 
standards  the  majority  of  the  time  while  iron  only  would  be  expected  to 
exceed  standards  occasionally  in  the  headwaters  of  the  reservoir. 

35-  During  periods  of  algal  blooms,  water  supply  treatment  costs 
could  be  increased.  The  magnitude  of  the  1973  algal  bloom  could  diminish 
the  recreational  potential  of  the  project.  The  other  simulated  algae 
blooms  would  probably  not  interfere  with  recreational  purposes  of  the 
lake  to  any  extent  greater  than  that  occurring  in  other  nutrient-rich 
Oklahoma  reservoirs  presently  receiving  heavy  recreational  use. 

36.  The  principal  study  indicated  pesticides  and  other  heavy 
metals  were  not  anticipated  to  vitiate  project  purposes  of  flood  "ontrol, 
municipal  and  domestic  water  supply,  and  general  recreation.  While 


13 


results  from  this  study  do  not  alter  this  conclusion,  it  is  important 
to  realize  this  conclusion  is  based,  in  part,  on  EPA’s  restricted  use 
of  potential  problem  pesticides  and  their  adsorption  and  sedimentation 
with  particulate  matter  as  the  suspended  solids  settle  out  in  the  head- 
water area  of  the  impoundment . Although  turbidity  and  suspended  solids 
concentrations  are  expected  to  be  greater  in  the  water  column  at  the 
lower  pool  elevation,  sampling  studies  on  Lake  Thunderbird  and  Lake 
Eufala  in  the  previous  study  indicated  pesticide  and  heavy  metal  concen- 
trations in  the  water  column  decreased  down  the  length  of  the  reservoir. 
The  sedimentation  of  these  constituents  may  occur  farther  into  the 
reservoir,  but  based  on  information  presently  available  pesticides  and 
heavy  metals  are  not  anticipated  to  be  water-quality  problems  for 
established  project  purposes. 

37-  Total  dissolved  solids  are  expected  to  exceed  water-quality 
standards  throughout  the  proposed  Arcadia  Lake.  Standards  for  this 
constituent  are  based  on  aesthetic  rather  than  toxicological  effects. 
Fecal  coliform  standards  for  primary  contact  recreation  were  not  vio- 
lated in  the  simulated  area  near  the  dam  in  either  the  1970  or  1973 
computations.  Standards,  however,  might  be  violated  in  the  upper  third 
of  the  reservoir  during  most  storm  events  and  throughout  the  impoundment 
during  major  storm  events. 

38.  Specific  conclusions  concerning  the  water  quality  of  the 
proposed  Arcadia  Lake  at  pool  el  305  m msl  are  based  on  information 
obtained  from  several  approaches  and  a state-of-the-art  evaluation  of 
these  data.  While  the  use  of  several  approaches  provides  corroborating 
evidence  for  each  specific  conclusion,  these  conclusions  are  not  un- 
equivocal. There  is  a degree  of  uncertainty  associated  with  each  ap- 
proach and  varying  degrees  of  uncertainty  associated  with  each  parameter 
being  evaluated.  For  example,  simulations  of  the  thermal  regime  expected 
in  the  proposed  Arcadia  Lake  are  based  on  physical  processes  that  are 
reasonably  well  understood.  Simulations  of  the  biological  processes, 
however,  are  based  on  empirical  descriptions  that  are  not  as  well  un- 
derstood. With  the  use  of  loading  analyses,  algal  bioassays,  and 
analysis  of  data  from  surrounding  impoundments,  conclusions  on  the 
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eutrophication  potential  can  be  made  with  a relatively  high  degree  of 
certainty.  Based  on  the  hydraulic  residence  time,  sensitivity  analyses, 
and  decay  processes,  conclusions  on  coliform  bacteria  are  also  rela- 
tively certain  for  the  area  near  the  dam.  This  is  not  true  for  the 
headwater  area  of  the  proposed  impoundment.  Sediment  transport  phenomena 
are  poorly  understood  and  mathematical  simulations  of  sediment  transport 
provide  only  gross  estimates  of  expected  conditions.  Particle-size 
distribution,  density  differences  between  inflowing  stream  and  impound- 
ment, inflow  entrainment,  wind  mixing,  local  morphometric  characteristics, 
and  other  factors  all  have  an  influence  on  the  rate  and  location  of 
sedimentation.  Since  coliforms  are  known  to  adhere  to  sediment  particles, 
their  distribution  will,  in  part,  be  dictated  by  the  sediment  distribu- 
tion. As  evidenced  by  delta  formation  in  the  headwater  areas  of  most 
impoundments , much  of "the  sedimentation  occurs  in  the  headwater  area; 
therefore,  most  of  the  coliform  colonies  will  also  be  deposited  there. 

The  coliform  colonies  that  remain  attached  to  suspended  particles  die  at 
an  exponential  rate. 

39-  While  additional  pesticide  and  heavy  metal  data  were  not  col- 
lected during  this  study,  the  lower  pool  elevation  may  modify  some  of 
the  conclusions  on  contaminants  stated  in  the  principal  study.  Because 
of  these  potential  modifications,  the  pesticide  and  heavy  metal  data 
were  reevaluated  in  light  of  conditions  expected  at  the  lower  pool  ele- 
vation. The  majority  of  pesticides  and  heavy  metals  entering  Arcadia 
Lake  are  expected  to  be  adsorbed  on  suspended  particles.  As  with  coli- 
forms, in-lake  pesticide  and  heavy  metal  concentrations  are  expected  to 
be  significantly  reduced  by  sedimentation  in  the  headwater  areas  of  the 
impoundment.  Unlike  coliforms,  however,  some  pesticides  and  heavy  metals 
may  undergo  complex  chemical  transformations  that  may  increase  the  solu- 
bility of  the  constituent  under  certain  conditions.  These  transforma- 
tions are  affected  by  a variety  of  conditions  such  as  pH,  oxygen  regime, 
redox  potential,  chelating  agents,  and  other  factors  that  are  difficult 
to  predict.  The  degree  of  certainty  of  the  conclusions  on  pesticides 
and  heavy  metals,  therefore,  is  reduced.  This  reflects  both  the  vari- 
ance in  the  data  base  and  the  current  state-of-knowledge  surrounding 
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the  environmental  chemistry  of  pesticides  and  heavy  metals.  Pesticides 
do  degrade  to  innocuous  compounds  hut  at  various  rates.  In  general, 
the  organophosphorus  compounds  degrade  more  rapidly  than  the  chlorinated 
hydrocarbons.  Partially  for  this  reason,  EPA  has  placed  restrictions  or 
banned  the  use  of  various  chlorinated  hydrocarbons.  Based  on  these 
restrictions,  the  concentrations  of  the  pesticides  indicated  as  potential 
problems  in  this  report  and  the  principal  study  report  are  expected  to 
diminish.  This  is  supported  by  their  reduction  in  the  water  columns  of 
Lake  Eufala  and  Lake  Thunderbird  from  the  headwaters  toward  the  dam  and 
by  a general  decreasing  concentration  trend  of  several  of  the  pesticides 
in  the  Deep  Fork  River  since  use  restrictions  were  imposed.  The  fish 
collected  in  the  Deep  Fork  River  clearly  violated  the  FDA  body  burden 
limit  for  -heptachlor  -epoxide  on  that  particular  sampling  date;  Since' 
aquatic  organisms  do  depurate  contaminants  and  analytical  anomalies  do 
occur  in  pesticide  determinations,  additional  samples  will  be  required 
to  determine  if  this  violation  was  the  norm,  an  analytical  error,  or  an 
improbable  occurrence.  Heptachlor  epoxide  was  not  detected  in  the  water 
column  or  the  sediments  of  the  Deep  Fork  River  near  Arcadia  during  the 
principal  study.  However,  heptachlor  was  found  in  the  water  but  not  the 
sediments  on  22  May  197I+ . Both  heptachlor  and  heptachlor  epoxide  were 
detected  in  the  Deep  Fork  River  during  the  1976  ACOG  sampling  program, 
but  concentrations  of  both  were  below  the  recommended  EPA  criteria  for 
the  protection  of  aquatic  life.^  While  various  pesticides  can  be 
indicated  using  gas  chromatography,  definitive  identification  requires 
further  analysis  using  mass  spectrometry.  It  is  possible  that  the  com- 
pound exceeding  the  FDA  criteria  for  fish  body  burden  is  not  heptachlor 
epoxide  but  some  compound  with  a similar  chromatographic  peak  residence 
time.  An  analysis  of  fish  collected  in  Oklahoma  reservoirs  as  part  of 
the  National  Pesticide  Survey  indicated  excessive  bioaccumulation  of 
pesticides  was  not  occurring  in  existing  impoundments  including  Lake 
Eufala  which  now  receives  the  runoff  from  the  Deep  Fork  River. 

1+0.  Mercury  can  exist  in  various  chemical  states,  has  known 
toxicological  characteristics,  and  has  been  measured  in  the  Arcadia 
watershed  and  Deep  Fork  River.  Mercury  concentrations  in  the  Deep  Fork 
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River  have  been  found  to  exceed  the  criterion  for  freshwater  aquatic 
life  but  do  not  exceed  the  criterion  for  drinking  water.  Unfortunately, 
the  criterion  for  freshwater  aquatic  life  is  below  the  detection  of  the 
analytical  equipment  used  by  the  USGS,  Oklahoma  Public  Health  Department, 
Oklahoma  State  University,  and  others  who  have  collected  water-quality 
data  on  the  Deep  Fork  River.  The  sensitivity  of  the  analytical  equip- 
ment used  by  the  Oklahoma  Public  Health  Department  is  properly  estab- 
lished to  detect  mercury  concentrations  for  drinking  water  criteria. 

The  criterion  for  freshwater  aquatic  life  is  nearly  two  orders  of  magni- 
tude lower  than  that  for  drinking  water  and  below  the  detection  limit  of 
the  analytical  equipment  used  for  analyses.  The  fish  collected  from  the 
Deep  Fork  River  did  not  exceed  the  FDA  body  burden  limit  for  mercury. 
However,  the  same  argument  applies  for  mercury  as  it  did  for  the  pesti- 
cides.. This  .fish .sample  may  represent  the  norm  or  it'  could  represent 
an  improbable  occurrence.  That  is,  additional  samples  might  indicate 
the  body  burden  limit  for  mercury  in  fish  was  exceeded.  Mercury  concen- 
trations in  storm  runoff  from  the  Deep  Fork  Basin  above  Belle  Isle  Lake 
have  been  found  to  exceed  both  drinking  water  and  freshwater  aquatic 
life  criteria. ^ Mercury  data  obtained  for  storm  runoff  and  under  ele- 
vated flow  conditions  in  the  Deep  Fork  River  near  Arcadia  by  the  USGS 
and  the  Oklahoma  State  University  probably  are  more  indicative  of  concen- 
trations to  be  expected  in  storm  event  loadings  to  Arcadia  Lake.  However, 
the  storm  runoff  data  from  above  Belle  Isle  Lake  indicate  that  the  ex- 
treme upper  Deep  Fork  Basin  may  be  a potentially  significant  source  of 
mercury. 

Conclusions  and  Recommendations 


Ul.  Based  upon  the  above  constraints,  the  data  base  existing  for 
this  study,  and  the  current  state  of  the  art  for  water-quality  analyses, 
specific  conclusions  concerning  the  water  quality  of  the  proposed 
Arcadia  Lake  are  as  follows: 

a.  The  proposed  Arcadia  Lake  would  be  thermally  stratified 
intermittently  during  the  summer  months.  No  permanent 
or  stable  temperature  stratification  is  expected  at  pool 
el  305  m. 
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b.  During  the  summer  and  early  fall  months,  release  tempera- 
tures may  be  about  T°C  warmer  than  natural  downstream 
temperatures . 

c_.  In  general.,  dissolved  oxygen  concentrations  will  be 

greater  than  5 mg/Jl  with  some  anoxic  periods  during  brief 
periods  of  intermittent  stratification. 

<1.  The  aerobic  environment  should  reduce  the  level  of  iron 
and  manganese  concentrations  near  the  dam.  However, 
manganese  is  expected  to  violate  drinking  water  standards 
the  majority  of  the  time  while  iron  may  occasionally  ex- 
ceed standards  in  the  headwater  area  of  the  impoundment. 
Total  dissolved  solids  are  expected  to  exceed  standards 
throughout  the  impoundment.  The  standards  for  iron, 
manganese,  and  TDS  are  based  on  aesthetic  rather  than 
toxicological  considerations  and  will  not  result  in  major 
water  supply  problems  if  considered  in  treatment  plant 
design. 

e_.  Pesticides  and  heavy  metals  are  not  anticipated  to  vitiate 
project  purposes.  Fecal  coliform  standards  for  primary 
contact  recreation  were  not  violated  near  the  dam  at  any 
time  during  the  simulations.  Coliforms  are  expected  to 
exceed  standards  in  the  headwater  area  of  the  impoundment 
and  may  be  high  in  the  upper  third  of  the  impoundment 
following  most  storm  events.  The  standards  may  be  ex- 
ceeded throughout  the  reservoir  following  major  storm 
events . 

_f . Two  or  three  major  algae  blooms  may  be  expected  each  year 
depending  upon  hydrometeorological  conditions.  It  is  ex- 
pected these  blooms  will  be  similar  in  magnitude  to  those 
occurring  in  Lakes  Keystone,  Thunderbird,  and  Carl  Black- 
well  and  may  cause  taste  and  odor  problems  in  the  water 
supply.  These  potential  problems  can  be  provided  for  in 
treatment  plant  design. 

Water-quality  requirements  for  project  purposes  of  flood 
control,  municipal  and  industrial  water  supply,  and 
general  recreation  are  expected  to  be  met  by  the  predic- 
ted water-quality  for  tne  proposed  Arcadia  Lake. 

1+2.  The  following  recommendations  are  made: 

a.  The  sampling  program  for  coliforms,  pesticides,  and  heavy 
metals  - particularly  mercury  - should  be  maintained 
throughout  the  planning,  construction,  and  operation 
phases  of  the  project.  Mercury  analyses  should  be  per- 
formed by  equipment  with  a detection  limit  well  below 
EPA  criteria  for  aquatic  life. 

b.  Prior  to  impoundment  additional  samples  of  fish  should  be 
collected  from  the  Deep  Fork  River  near  the  proposed 
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dams it e and  analyzed  for  body  burdens  of  mercury  and  the 
pesticides,  heptachlor  and  heptachlor  epoxide.  If 
heptachlor  or  heptachlor  epoxide  is  indicated  by  gas 
chromatography  to  be  in  significant  concentrations,  its 
identity  should  be  verified  by  mass  spectrometry.  If 
the  results  are  confirmed,  a more  thorough  evaluation  of 
its  source  and  potential  fate  and  consequences  in  the 
reservoir  should  be  made  prior  to  construction.  While 
present  evidence  does  not  justify  the  delay  of  advanced 
engineering  and  design,  the  CE  shcHd  be  prepared  to 
delay  plans  for  construction  if  the  results  of  the 
recommended  sampling  program  indicate  a potentially  sig- 
nificant problem. 

c_.  Recreation  sites  should  be  located  in  the  lower  half  of 
the  impoundment  to  minimize  potential  coliform  problems 
with  body  contact  recreation.  The  lower  half  represents 
a conservative  value  since  most  storm  events  are  expected 
to  impact  only  the  upper  third  of  the  proposed  impound- 
ment with  respect  to  coliforms-  violating  body  contact 
recreation  standards. 

cl.  The  selective  withdrawal  capability  should  be  maintained 
in  the  structure  to  permit  the  withdrawal  of  water  supply 
from  three  depths  within  the  impoundments.  These  three 
ports  should  probably  be  located  to  permit  the  withdrawal 
of  epilimnetic  water  (center  line  12  m above  the  bottom), 
metalimnetic  water  (center  line  8 m above  the  bottom),  and 
hypolimnetic  water  (center  line  1+  m above  the  bottom). 

The  hypolimnetic  port  is  in  addition  to  the  floodgates. 
During  algae  blooms,  water  may  be  withdrawn  from  the 
bottom  to  ameliorate  taste  and  odor  problems.  During 
stratification  periods,  surface  or  metalimnetic  waters 
may  be  withdrawn  to  minimize  manganese  and  iron 
concentrations . 
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Figure  3.  Release  temperatures  for  1970  at  the  lower  pool  elevation 
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Figure  k.  Release  temperatures  for  1973  at  the  lower  pool  elevation 
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Figure  J.  Algae  concentrations  predicted  during  1970  at  the  lower  pool  elevation 


Figure  8.  Algae  concentrations  predicted  for  1973  at  the  lower  pool  elevation 


Figure  9*  DO  isopleths  for  light  half-saturation  coefficient  during  1970 


Figure  11.  DO  isopleths  for  self-shading  coefficient  during  1970 


opleths  for  self-shading  coefficient  during  1973 


Figure  13.  Effect  of  light  half-saturation  coefficient  on  1970  algae  concentrations 


Figure  15.  Effects  of  self-shading  coefficient  on  algae  concentrations  during  1970 


Figure  17.  DO  isopleths  for  Secchi  disk  at  1.0  m during  1970 


sopleths  for  Secchi  disk  at  1.0  m during  1973 


Figure  19*  DO  isopleths  for  Secchi  disk  at  3.0  m during  1970 


Figure  20.  DO  isopleths  for  Secchi  disk  at  3.0  m during  1973 


21.  Effect  of  varying  light  penetration  on  1970  algae  concentrations 


Figure  22.  Effect  of  varying  light  penetration  on  1973  algae  concentrations 


Figure  25-  DO  isopleths  during  1970  at  pool  el  308 


Figure  26.  DO  isopleths  during  1973  at  pool  el  308 
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Figure  27.  1970  algae  concentrations  at  pool  el  308 


In  accordance  with  EE  70-2-3,  paragraph  6c(l)(b), 
dated  15  February  1973,  a facsimile  catalog  card 
In  Library  of  Congress  format  Is  reproduced  below. 
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